function [ theta ] = read_theta( n_mk, alpha )

    theta = zeros(size(n_mk,1),size(n_mk,2));
    for m = 1:size(n_mk,1)
        for k = 1:size(n_mk,2)
            top = n_mk(m,k)+alpha(k); 
            da = 0;
            for j = 1:size(n_mk,2)
                da = da + n_mk(m,j)+alpha(j);
            end
            theta(m,k) = top/da;
        end
    end

end

